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SUPER DECOUPLED LOADFLOW METODS 
ABSTRACT 

The complete patent specification involves the invention of 
six super decoupled loadflow methods. The Fast Super Decoupled 
Loadflow (FSDIi) method and two variants Transformed Fast 
Decoupled Loadflow (TFDLXB and TFDLBX) methods are known ones and 
included are their novel versions (NFSDL, NTEDLXB, NTFDLBX) . 
These are the best versions of many simple variants with almost 
similar performance. These methods specifically involved the use 
of the following invented technics in the prior art method. 

1. Gain matrices of all the six methods are different and they 
can be determined as described in the specification. 
Obviously the gain matrices of Pe-subproblem can be defined 
unsymmetrlcal . 

2. In all the six methods, rotation angles are restricted to 
the maximum value of -36 degrees from nonlinearity 
considerations. 

3. Slack-start for any decoupled loadflow method for efficiency 

4. Modification of real power mismatches at PV-nodes according 
to relations (12) and (15) of the specifications for the 
methods FSDL, TFDLXB, TFDLBX. 

5. ' Computation of angle correction and updating for PV-nodes 

along with the voltage magnitude corrections at PQ-nodes for 
NFSDL, NTFDLXB and NTFDLBX methods. In the back-substitution 
part of the solution of this subproblem involving angle 
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corrections at PV-nodes and voltage magnitude corrections at 
PQ-nodes, skip factor elements corresponding to rows and 
columns of PV-nodes in the calculation of voltage magnitude 
corrections . 

The specification deals only with unadjusted super decoupled 
loadflow solution methods and does not describe any applications. 
However the computer algorithms can be appropriately modified 
for adjustments and/or applications such as state estimation, 
contingency analysis etc. 

Also described and claimed is the invention of two compact 
storage schemes for gain matrices ofthe FSDL method and two 
efficient procedures for node type switching implementations in 
the FSDIi method. 

Dated this 7 th day of November 1993 




Signature of the inventor - s.B.Patel 
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THE CANADIAN PATENT ACT 



THE COMPLETE SPECIFICATION 



SUPER DECOUPLED LOADFIiOW METHODS 



Sureshchandra Bhagubhai Patel, Indian Nat:ional, 
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The following specification particularly describes and ascertains 
the nature of this invention and the manner in which it is to be 
per f ozrmed . 
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This inven-tion rela-tes -bo 'the super decoupled me-thods o£ 
s'teady-s'ta'te loadflow analysis of power sys-tem, and compu-ter 
algorit:hins for carrying out -these methods. 

A computer algorithm processes raw information to yield 
useful information. A chemical process processes raw material to 
yield useful material. Useful information is well recognized as 
technological product in the modern age of infomnation 
technology. Therefore, the computer algorithm like chemical 
process yields useful product and it is the useful art. 

INTRODUCTION 

Iioadflow in power system studies is the most basic frequen- 
tly performed steady state analysis of an electrical power net- 
work. The loadflows are performed in system planning, operational 
planning, and operation control. Of various methods proposed over 
the past four decades, Stott's Fast Decoupled Loadflow (FDI.) [1] 
has gained wide acceptance both for off- and on- line applica- 
tions. However, it is known to suffer poor convergence for 
systems having high R/X ratio branches (discussions in [1] & 
[2]) . 

Dy Liacco and Ramarao (discussion of [2]) as well as Deckmann 
et. al. [3] have developed circuit transformation schemes to 
avoid such difficulties. The two schemes however do not provide 
consistent improvement in convergence [ 4 ] . Also they are not 
general and are suitable where the system has a small number of 
troublesome R/X ratio branches. 

A recently introduced super decoupled approach [5,6,7] is 

more general and reliable. Rotation operators applied to the 

complex node injections and the corresponding admittance values 

2 
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that relate the above to the system state variables, transform 
the network equations such that branch admittances appear to be 
almost entirely reactive. Thus, better decoupling is realized. 
However, the super decoupled algorithms are not efficient and are 
not very convenient in contingency analysis. A recent modifica- 
tion to the FDL method for networks with high R/X ratio branches 
[4] is in fact variation of the super decoupled approach with a 
weakness that it takes higher number of iterations for normal 
cases . 

General-purpose version of the Fast Decoupled Loadflow 
(GFDL) method [8] proposes simple modifications to the classical 
FDL method. The convergence is much improved in the presence of 
large R/X ratio lines. The two versions of FDL method are put in 
a better theoretical framework by Monticelli et. al.[9]. The GFDL 
method is nothing but an experimental investigation of the 
original observation made about the behaviour of the FDL method 
by this inventor as back as 1985 in an unaccepted research paper 
communicated to IEEE (New York) . 

The critically coupled loadflow methods are based on the use 
of the overlap update rule [10]. The disadvantage of the methods 
is that they involve the solution of a larger number of equa- 
tions at each iteration. Moreover the tests conducted for the 
paper CIO] avoided the increase of R/X ratio of branches connected 
to the FV-nodes. Therefore, conclusions of the paper are not 
based on exhautive testing. Test show that the invented six 
versions are the best loadflow methods. However, simply modified 
versions including hybrids of the six invented methods can also 
have closely similar performance for any given system. 



3 
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The laliest: Novel Decoupled liOadflow Met:hod[ll] can simply be 
derived from published [ 12 ] and lately claimed invetion of this 
patent specification. In fact the idea involved in the method of 
[11] was revealed by this inventor in extensive communication for 
about four years with Institution of Electrical Engineers of 
England. Also the idea was investigated along with the modifica- 
tion of the real power mismatches at PV-nodes and only the best 
method was reported in reference [12]. Moreover the unsymmetrical 
gain matrix definition for the Pesubproblem is obvious from the 
definition of the factor (K) used to modify real power mismatches 
of the PV-nodes [ 12 ] . However the method of [11] is not reliable 
for systems having PV-nodes and it has no flexibility of tuning 
rotation angle as was originally known to this inventor. Possibly 
this may be the reason that the author of [11] did not comapre 
his method with the Fast Super Decoupled Loadf low [ 12 ] and 
Transformation based Fast Decoupled liOadf low[ 13 ] methods. 

Keywords : Power systems, Loadflow, Matrices, solution of 

simultaneous equations. 

The invention will now be described, wherein the following 
symbols are used : 

Y =G +jB : (p-q)-th element of nodal admittance 
pq pq pq matrix formed excluding shunts 

y = g +jb : total shunt admittance at node p 
p p p 

V =»e +jf ="V : complex voltage at node p 
P P P P P 

<\® ' ' voltage angle, magnitude corrections 

P P 
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Ae , : real, imeginary coroponen-t corrections of voltage 
P P 

P + jQ : net nodal injected power calculated 
P P 

+j /5iQ s nodal power residue (mismatch) 
P P 

RP +jRQ : modified nodal power residue 
P P 

m : number of PQ-nodes 

k : number of PV-nodes 

n=m+k+l : total number of system nodes 

q>.p : q is the node adjacent to node p excluding the case q 

[ ] : indicates enclosed variables to be vector or matrix 



SUPER DECOUPLED LOADFLOW (SDL) METHOD { The prior art) 

There are two versions of the Super Decoupled Loadflow 
method. These are based on the XB- and BX- versions of the Fast 
Decoupled Loadflow (FDL) method. The one based on the XB-version 
was developed in 1985 [7] and the other based on the BX-version 
is obvious from the BX-version of the FDL method developed in 
1989 [8]. These methods involve the iterative solution of system 
of equations (1) and (2) . Unlike FDL model the gain matrices in 
its, transformed version are unsymmetrical because mostly 
different rotations are required to be applied at the terminal 
nodes of a branch. Symmetrical gain matrices can be obtained by 
the Haley and Ayres technique [7] of applying average of the 
rotations at the terminal nodes of a branch to the branch 
admittance (appendix) . 
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SDL based on t:he XB-verslon of the FDIj (SDLXB) 

The SDLXB-version is described by the following equations 
(1) to (10) . 

[RP] = [Y'] [ £^0] (1) 
[RQ] = [Y"] [ ti^V] (2) 

Wherein each element of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of [RP] at pv-nodes 
is given by (5) . 

RP = ( liP Cosp + Sin/a )/V (3) 
P P P P P P 

RQ = (-£iP Sinp + jS,Q Cosp )/V (4) 
P P P P P P 

RP = t^P /V (5) 
P P P 

Trigonometric functions in (3) and (4) and elements of [Y'] and 
[Y"] are given by (6), (7) and, (8), (9) and (10). 



P PP V PP PP 



P 



Cosp = -B / /g + (6) 



r~2 2 

Sinp = -G / /G + B (7) 
P PP V PP PP 

Y' = -1 / X and Y" = -B (8) 

pg pq pq pq[ 

Y' =JE-^' ^^'^ ^" = -2b'+ X (9) 

PP f-*f pq pq p t^r pq 

b' = b Cosp or = b (10) 
P P P P 

Where X is the transformed branch reactance defined in appendix 
pq ^ 

by eqatlon (41) and B is the corresponding transformed element 

pq. 

of the susceptance matrix. 

6 
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SDIi based on the BX-version of the FDL (SDLBX) 

The SDLBX-version differs from the SDLXB-version only in 
relation (8) as given by relation (11) . The SDIiBX-version 
consists of relations (1) to (7) , (11) , (9) and (10) . 

Y' = -B and Y" = -1 / X (11) 

pq pq pq pq 

INVENTED SUPER DECOUPLED LOADFLOW METHODS 

Six invented super decoupled loadflow models are described 
in this section. Fast Super Decoupled and two versions of ^ the 
Transformed Fast Decoupled models are Known ones [ 12 , 13 ] and 
included are their novel versions. 

Fast Super Decuopled Iioadflow (FSDL) 

The FSDL method involves the iterative solution of the 
system of equations (1) and (2) . The model cosists of relations 
(1) to (4), (6), (7), (12) to (16), (9) and (10). 

RP - tJ? /(K V ) (12) 
P P P P 

Elements of [Y'] and [Y"]. and the multiplier K in (12) are given 

P 

by : 

Y ' = -Y for branch r/x ratio < 2.0 

pq pq 

= -{B +0.9(Y -B )> for branch r/x ratio > 2.0 
pq pq pq 

= -3pq for branches connected between 

two PV-nodes or a PV-node and 
the slack-node (13) 

Y" = -Y for branch r/x ratio < 2.0 

pq pq 

= -{B +0.9(Y -B )} for branch r/x ratio > 2.0 (14) 
pq pq pq 

7 
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K = Al3Solut:e (B /Y' ) (15) 
P PP PP 

KKp = Absolutie (Y'pp/Y"pp) (16) 



Branch adini1:t:ance inagnit:ude In (13) and (14) is of the same 
algebraic sign as it:s susceptiance . Elementis of the two gain 
matrices differ in that diagonal elements of [Y"] additionally 
contain the b' values given by equation (10) and in respect of 
elements corresponding to branches connected between two PV-nodes 
or a PV-node and the slack-node. The factor 0.9 in the relations 
(13) and (14) is tuned only for rotation limited to -3 6 degrees. 
With different Limiting Rotation Angle (UlA), it needs to be tuned 
again. In two simple variations of the FSDL method, one is to 
make Y" = Y' and the other is to make Y' = y" 

pq pq pq pq 

Transformation based Fast Decoupled Ijoadflow (TFDL) 

The TFDL model is similar to the FSDL model. They differ 
only in the definition of the gain matrices. 
The TFDL(XB-version) (TFDLXB) Loadflow 

The TFDLXB-version consists of relations (1) to (4), (12), 
(15), (6), (7), (17), (18), (19), (9) and (10). 

Y'pq =f-Bpq for branches connected between two PV- 

nodes or a PV-node and the slack-node 

-1 / X for all other branches (17) 

pq 

Y" = -B (18) 
pq pq 

KK = Absolute (Y' / ( X V )) (19) 

p PP -f^7 pq 

8 
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The TFDL(BX-version) (TPDLBX) Loadflow 

The TFDLBX-version consists of relations (1) to (4), (12), 
(15), (6), (7), (20), (21), (22), (9) and (10). 

Y'pq =p^pq for branches connected between two PV- 

nodes or a PV-node and the slack-node 

-B for all other branches (20) 

pq 

Y" = -1 / X (21) 

pq pq 

KK = Absolute (Y' /( )) (22) 

p pp i-*!* pq 

Where X is the transformed branch reactance defined in appendix 

pq 

by eqation (41) and B is the corresponding transformed element 

pq 

of the susceptance matrix. 



The factor KK is to be multiplied to the real power mis- 
P 

match at a PV-node switched to PQ-type in node type switching 

implementation. From general considerations, K and KK are 

P P 

restricted to the minimum value of 0.75 for FSDL, TFDI<XB and 
TFDLBX methods. However they can be tuned for the best possible 
convergence for any given system. 



In all the SDLXB, SDLBX, FSDL, TFDLXB and TFDLBX models 
[Y'] and [Y"] are real, sparse, symmetrical and built only from 
network elements. Since they are constant, they need to be 
factorized once only at the start of the solution. Equations (1) 
and (2) are to be solved repeatedly by forward and backward 

substitutions. 

. • 9 
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[Y'] and [Y"] are of the same dimensions (m+k)*(m+k) when 
only a row/column of the slack-node is excluded and both are 
triangularized using the same ordering regardless of the node- 
types. For a row/ column corresponding to a PV-node excluded in 
[Y"], use a large diagonal to mask out the effects of the off- 
diagonal terms. When the node is switched to the PQ-state the 
row/column is reactivated by removing the large diagonal. This 
technique is especially useful in the treatment of PV-nodes in 
the matrix [Y"]. 

It is invented to make this technique efficient while 
solving (2) for by skipping all PV-nodes and factor elements 
with indices corresponding to PV-nodes. In other words efficiency 
can be realized by skipping operations on rows/columns 
corresponding to PV-nodes in the forward-backward solution of (2) 
for J^V. This has been implemented and the time saving of about 
4% of the total solution time (including input/ output) could be 
realized in 14-14 iterations required to solve 118-node system 
with the uniform R-scale factor 4 applied. The time saving has 
been assessed on PC-XT. It should be noted that the same indexing 
and addressing information can be used for the storage of both 
the matrices as they are of the same dimensions and sparsity 
structure . 

Novel Fast Super Decoupled Loadflow (NFSDL) 

I 

The NFSDL method involves the iterative solution of the 
following system of equations (23) and (24) . 



[RP] = [Yl] [ £^0] 



(23) 




(24) 
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Wherein each element: of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of [RP] at PV-nodes 
is given by (5). In (24) PV-nodes are assumed to be grouped and 
numbered after all PQ-nodes. 

Elements of [Yl] and [Y2] are given by : 



pq 



-Y 

pq 

-{B +0.9(Y -B 

pq pq pq 

-Bpq 



for branch r/x ratio ^ 2.0 
) } for branch r/x ratio > 2.0 



for branches connected between 
two PV-nodes or a PV-node and 
the slack-node (25) 



Y2pq = 



Gpq for branches connected to a PQ-node 

and a PV-node 



-Ylpq for branches connecting two PQ-nodes 



Yl 



pq 



for branches connecting two PV— nodes 



(26) 



Yl 

PP 
Y2pp 



pq 



^b^>-Ylpp' 



-B 



PP 



for PQ-nodes 



for PV-nodes 



(27) 



(28) 



The NFSOIi model consists of relations (23), (24), (3) to 
(7), (25) to (28), and (10). Branch admittance magnitude in 
(25) is of the same algebraic sign as its susceptance. 

Novel Transformed Fast Decoupled Iioadflow (NTFDIi) 



The NTFDL model is similar to the NFSDL model. They differ 

only in the defination of gain matrices. 

11 
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The NTFDL(XB-version) (NTFDLXB) Loadflow 

The NTFDIiXB model consis'bs of relatiions (23), (24), (3) -to 
(7), (29) to (32), and (10). 



Ylpq = 



Y2pq 



Yl 



PP 



-Bpq 
-1 / X 

I — pq 



= rGpq 



pq 



-B 



pq 



*T-»-p pq 



for branches connecked between two PV- 
nodes or a PV-node and the slack-node 

for all other branches (29) 



for branches connected to a PQ-node 
and a PV-node 



for branches connecting two PQ-nodes 



for branches connecting two PV-nodes (30) 



(31) 



Y2 



PP 



2b'+B 
P PP 



PP 



for PQ-nodes 



for PV-nodes 



(32) 



Where X is the transformed branch reactance defined in appendix 
pq 

by eqation (41) and B is the corresponding transformed element 

pq 

of the susceptance matrxx. 

The NTFDL(BX-version) (NTFDLBX) Loadflow 

The NTFDIiBX model consists of relations (23), (24), (3) to 
(7), (33) to (36), and (10). 



Ylpq = 



-Bpq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



-B 



for all other branches 



(33) 



pq 



12 
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Y2pq 



Gpq 



1/ X 



pq 



for branches connected to a PQ— node 
and a PV-node 

for branches connecting two PQ-nodes 



. B for branches connecting two PV-nodes (34) 

L.pq 



Yl 



PP 



Y2 



PP 



= 21 
i-»-p pq 



2b' + 51-1/ X 

p t-*p pq 



-B 



(35) 



or = 2b '+1/5? for PQ-nodes 
P PP 

for PV-nodes (36) 



PP 



In the NFSDL NTFDIiXB and NTFDLBX methods [Yl] and [Y2] are 
real, sparse, symmetrical and built only from network elements. 
Since they are constant, they need to be factorized once onl^ at 
the start of the solution. Equations (23) and (24) are to be 
solved repeatedly by forward and backward substitutions. 

[Yl] and [Y2] are of the same dimensions (m+k) * (m+k) when 
only a row/column of the slack— node is excluded and both are 
triangularized using the same ordering regardless of the node- 
types. It should be noted that the same indexing and addressing 
information can be used for the storage of both the matrices as 
they are of the same dimensions and sparsity structure. Unlike 
[Y"], It is to be noted that all the PV-nodes are also active in 
[Y2]. Therefore, the novel methods would pose a difficult problem 
of node-type switching. 

ft is restricted to the maximum of -3 6 degrees from 
P 

nonlinearity considerations for FSDL, TFDXjXB, TFDLBX, NFSDIi, 
NTFDLXB and NTFDLBX methods. However it can be tuned for the 
best possible convergence for any given system and it can also be 



vertually unrstricted -90 degrees. 

13 
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It.era'tion Scheme 

The basic i'tera'tive scheme for solving "the proposed models 
is -to solve for [^9] "to update [6] and "then solve (2) or (24) . 
This is -the block Gauss-Seidel approach. The scheme is block- 
successive which imparts increased s-tabili-ty to the solution 
process. This in turn improves convergence and increases 
reliability. Cycling behaviour in the iterative process of the 
standared iteration scheme of [1] can be avoided by simply using 
different convergence tolerances for the real and reactive power 
mismatches as suggested in [9]. The cycling behaviour and its 
remedy as above were originally observed by this inventor as back 
as in 1985 in an unaccepted paper submitted to IEEE and prepared 
from research conducted at the university of Roorkee. However, 
algorithms are given for strictly successive iteration scheme of 
[8] . 

APPENDIX 

Transformation of Branch Admittance 

The branch admittance transforamation for symmetrical gain 
matrices of the TFDLXB, TFDLBX, NTFDIiXB and NTFDLBX methods is 
given by the following steps : 

1. Compute ; ft = arctan (G /B ) and 

pp pp pp 

(37) 

p = arctan (G /B ) 
qg gg qq 

2. Compute the average of rotations at the terminal nodes (p 
and q) of a branch : 

0 = (P + p ) / 2 (38) 
av pp qq 

14 
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3. Compare j3 with the Limiting Rotation Angle (LRA) and let 

av 

/3 to be the smaller of the two : 
av 

p = minimum (p , ISRA) (39) 
av av 

4. Compute transformed pq-th element of the admittance matrix : 

G +jB = (Cosp +jsinp ) (G +jB ) (40) 
pq pq av av pq pq 

5. Note that the transformed branch reactance is : 

2 2 

X = B / (G +B ) (41) 
pq pq pq pq 



and similarly 



2 2 

X = B / (G +B ) (42) 
PP PP PP PP 



AIjGOR1THM-1 (The prior art) 

Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (18, IV) Iterative Scheme : 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set specified 
voltage magnitudes at PV-nodes and 1.0 p.u. voltage 
magnitudes at PQ-nodes. Set all the node angles equal to 
that of the Slack-node angle. This is referred to as the 
flat-start [1]. 

(b) Initialize iteration counts ITRP=ITRQ=r=0 

(c) Form nodal admittance matrix 

15 
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(d) Form (in+k) by (m+k) size matrices [Y'] and [Y"] of (1) and 
(2) respectively each in a compact storage exploiting 
spars ity. 

(i) In case of SDLXB-method, the matrices are formed usig 
relations (8) , (9) and (10) . 

(ii) In case of SDLBX -method, the matrices are formed using 
relations (11), (9) and (10). 

(e) In [Y"] matrix, replace diagonal elements corresponding to PV 
-nodes by very large value (say 10. 0**10). In case [Y"] is of 
dimension (m by m),thls is not regired to be performed. 
Factorize [Y'] and [Y"] using the same ordering regardless 
of the node-types. In case [Y**] is of dimension (m by ra) , it 
is factorized using different ordering than that of [Y']. 

r r 

(f) Compute residues (PQ- and PV— nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( ^ ) , 
proceed to step (n) . Otheirwise follow the next step. 

(g) Compute the vector of modified residues [RP] using (3) for 
PQ-nodes and using (5) for PV— nodes. 

(h) Solve (1) for voltage angle corrections and update voltage 
angles using egn. (43) . 

(i) Increament the iteration count ITRP=ITRP+1 and 
r= (ITRP+ITRQ) /2 . 

r r 
(j) Compute residues /S^P (PQ- and PV-nodes) and f^Q (at PQ- 
nodes only) . If all are less than the tolerance ( ^ ) , 

proceed to step (n) . Otherwise follow the next step. 

16 
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(k) Compute the vector of modified residues [RQ] using (4) for 
PQ-nodes only. 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ-nodes only using eqn. (44) . 

(m) Increament the iteration count ITRQ=ITRQ+1 and 

r=(ITRP+ITRQ)/2 . Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

The SDIiXB and SDIiBX algorithms differ only in step-d defining 
gain matrices in the above solution steps. Fig.l is the flowchart 
of algorithm-1. 

AIjGORITHM-2 (containing invented steps) 

Fast Super Decoupled I.oadflow Algorithms (solution Steps) 
Strictly Successive (le, IV) Iterative Scheme | 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set all the 
nodes voltage magnitudes and angles equal to those of the 
Slack-node. This is referred to as the slack-start. 

(b) Initialize iteration counts ITRP=ITRQ=r=0 

(c) Form nodal admittance matrix 

(d) Form (m+k) by (m+k) size matrices [Y'] and [Y»«] of (1) and 
(2) respectively each in a compact storage exploiting 
spars ity- 

(i) In case of FSDL-method, the matrices are formed using 

relations (13), (14), (9) and (10). 

17 



2107388 



(ii) In case of TFDLXB-method , the matrices are formed uslg 
relations (17), (18), (9) and (10). 

(iii) In case of TFDLBX-method , the matrices are formed 
using relations (20), (21), (9) and (10). 

(e) In [Y"] matrix, replace diagonal elements corresponding to PV 
-nodes by very large value (say 10. 0**10). In case [Y"] is of 
dimension (m by m) , this is not reqired to be done. 
Factorize [Y' ] and [Y"] using the same ordering regardless 
of the node-types. In case [Y"] is of dimension (m by m) , it 
is factorized using different ordering than that of [Y']. 

r r 

(f) Compute residues (^P (PQ- and PV-nodes) and [^Q (at PQ- 

nodes only) . If all are less than the tolerance ( ) , 
proceed to step (n) . Otherwise follow the next step. 

(g) Compute the vector of modified residues [RP] using (3) for 

PQ-nodes and using (12) and (15) for PV-nodes. In (3) )3 is 

P 

restricted to the maximum of -3 6 degrees (or any other 
angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degree) . 

(h) Solve (1) for voltage angle corrections and update voltage 
angles using relation (43) . 

r (r-1) r 

e = e + Ke (43) 

P P P 

(i) Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 
iteration count ITRP=ITRP+1 and tc^ l^s:.'rp-P -^r X.TR^) ^ 2. ' 

IV 
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r r 
(j) Compute residues f\,P (PQ- and PV-nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( 6. ) , 

proceed to step (n) . Otherwise follow the next step. 

(k) Compute the vector of modified residues [RQ] using (4) . In 

(4) J2» is restricted to the maximum of -3 6 degrees (or , any 
P 

other angle) from nonlinearity considerations. It can be 
even unrestricted for any given system (-90 degree) . 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ-nodes only as given by relation 
(44) . While solving (2) , skip all the rows and columns 
corresponding to PV-nodes to increase the efficiency. 

V = V + CsiV (44) 

P P P 

(m) Increament the iteration count ITRQ=ITRQ+1 and 
r=(ITRP+ITRQ)/2 . Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical 
(le, IV) iteration scheme of reference [1]. This algorithm is for 
solving loadflow problem formulated in polar coordinates. The 
solution of the loadflow model formulated in rectangular 
coordinates also involve similar steps with minor modifications. 
Gain matrices of the invented methods are the same for both 
polar- and rectangular-coordinate formulations of the loadflow 
problem. Fig. 2 is the flowchart of algorithm-2 . The FSDL-method, 
TFDLXB-method and TFDLBX-method differ only in step-d defining 
gain matrices in the above algorithm-2. 
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AliGORITHM-3 (invclude invented steps) 

Novel Fast Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (1©, IV) Iterative Scheme | 

(a) Assign an initial approximate solution. If better solution 
estimate is not available, set all the nodes voltage 
magnitudes and angles equal to those of the Slack-node. This 
is referred to as the slack-start. 

(e) Initialize iteration counts ITRP=ITRQ=r=0 

(c) Form nodal admittance matrix 

(d) Form (m+k) by (m+k) size matrices [Yl] and [Y2] of (23) and 
(24) respectively each in a compact storage exploiting 
spars ity. 

(iv) In case of NFSDL-method , the matrices are formed using 
relations (25) to (28), and (10). 

(v) In case of NTFDLXB-method , the matrices are formed 
using relations (29) to (32) and (10) . 

(vi) In case of NTFDLBX-method , the matrices are formed 
using relations (33) to (36) and (10). 

(e) Factorize [Yl] and [Y2] using the same ordering regardless 
I of the node-types. 

r r 

(f) Compute residues fj^P (PQ- and PV-nodes) and |s,Q (^^ 

nodes only). If all are less than the tolerance ( ) , 
proceed to step (n) . Otherwise follow the next step. 

20 
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(g) Compute the vector of modified residues [RP] using (3) for 

PQ-nodes and using (5) for PV-nodes. In (3) p is restricted 

P 

to the maximum of -36 degrees (or any other angle) from 
nonlinear ity considerations. It can be even unrestricted for 
any given system (-90 degrees) . 

(h) Solve (23) for voltage angle corrections and update voltage 
angles using relation (43) . 

(i) Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 
iteration count ITRP=ITRP+1 and r= (ITRP+ITRQ) /2 . 

r r 
(j) Compute residues ^P (PQ- and PV-nodes) and t^Q (at PQ- 
nodes only) . If all are less than the tolerance ( €. ) , 
proceed to step (n) . Otherwise follow the next step. 

(k) Compute the vector of modified residues [RQ] using (4) for 

PQ-nodes. Compute [RP] for PV-nodes using (5). In (4) p is 

P 

restricted to the maximum value of -36 degrees (or any other 
angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degrees) . 

(1) Solve (24) for voltage magnitude corrections at PQ-nodes 
and voltage angle corrections at PV-nodes. Update voltage 
magnitudes of PQ-nodes using relation (44) and angles of PV- 
nodes using relation (43) . In the back substitution part of 
the solution of (24) , skip all the PV-nodes in the 
calculation of voltage-magnitude corrections of PQ-nodes. 
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(tn) 



Increament 



the 



iteration 



count 



ITRQ=ITRQ+1 



and 



r={ITRP+ITRQ)/2 . Proceed to step (f) 
(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical 
(le, IV) iteration scheme of reference [1]. This algorithm is for 
solving loadflow problem formulated in polar coordinates. The 
solution of the loadflow model formulated in rectangular 
coordinates also involve similar steps with minor modifications. 
Gain matrices of the invented methods are the same for both 
polar- and rectangular-coordinate formulations of the loadflow 
problem. The NFSDL-method , NTFDLXB-method and NTFDI.BX-method 
differ only in step-d defining gain matrices in the above 
algorithm-3. Fig. 3 is the flow-chart of algorithm-3. 

ALGORITHMS using GliOBAI. CORRECTIONS 

The algorithms-1, -2 and -3 of above involve increamental ( or 
local) corrections. All the above algorithms can be organised to 
produce corrections to the initial estimate solution. It involves 
storage of the vectors of modified residues and replacing the 
relations (3) , (4) , (5) , (12) , (43) and (44) respectivly by 
(45) , (46) , (47) , (48) , (49) and (50). Superscript 'O' in relations 
(49) and (50) indicates the initial solution estimate. 



r r r r 

RP = ( fc^P Cosfi + tiSi )/V + HP 

p P P P P P ] 



(r-1) 



(45) 



P 




(r-1) 



(46) 



P 
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r r r 

RP = j^P /V + RP (47) 

p p p p 

r r r (r-l) 

RP = IsP / (K *V ) + RP (48) 
P P P P P 

r Or 

e = e + (49) 
p p P 

r Or 
V = V + £n,V (50) 

P P P 



RECTANGULAR COORDINATE FORMULATIONS OF THE INVENTED LOADFLOW 
METHODS 



This involves following changes in the equations describing 
-the loadflow models formulated in polar coordinates. 



(i) Replace e and ^9 respectively by f and ^f in equations 
(1), (23), (24), (43) and (49). 



(ii) Replace V and ^s^V respectively by e and fv^e in equations 
(2), (24), (44) and (50). 



(iii) Replace V by e or e in equations (3), (4), (5), X12) , 

p p s 

(45) , (46) , (47) and (48) . The subscript 's' indicates the 
Slack-node variable. 



(iv) After calculation of corrections to the imeginary part of 

complex voltage ( of PV-nodes and updating the 

imaginary component (f) of PV-nodes, calculate real 
component by : 

e = - f^ (51) 

p p (specified) p 
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COMPACT STORAGE SCHEMES FOR GAIN MATRICES OF THE FSDL METHOD 

Two compact matrix storage schemes for gain matrices of the 
FSDL method are described in the following. 

(1) When network shunts are ignored (FK=0.0), [Y"] becomes the 
submatrix of [Y']. In this scheme, [Y'] is f actor ized using 
optimal ordering regardless of node-types along with ~ the 
following procedure : 

1. Identify PQ— nodes in the composite path of all PV-nodes. 

2 . Reload and ref actor all row/columns of these PQ-nodes using 
partial refactorization method PRl of [13]. Store these 
factors separately . 

3. When solving for jj^V using factors of [Y'] : 

(i) all PV-nodes and factor elements with indices 
corresponding to PV-nodes are skipped. Zeroing of reactive 
power mismatches at PV-nodes is not required. 

(ii) for PQ-nodes of the path, separately stored factors as 
in step-2 above are used along with the other factors of 
[Y'] in forward-backward operations. 

The above procedure is implemented. For 118-node system, 
there are only 18 PQ-nodes in the composite path of all 53 PV- 
nodes. Each of the gain matrices [Y'] and [Y"] is having 117- 
row/ columns for 118-node system, and the total of 2 34 -row/ columns 
are to be f actorized and stored for conventional 2-matrix 
solution. Threfore the proposed second scheme can be said to^ be 
saving 42% (99/234=0.42) of the factorization effort and core 
memory. This saving rises to the maximum 50% for systems without 
PV-nodes . 

24 
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(2) When network shunts are ignored (FK=0.0), [Y"] becomes the 
submatrix of [Y']. In the single matrix version of the FSDL 
method only [Y'] is factorized and stored. This is achieved 
exactly by placing PV-nodes after all PQ-nodes and by allowing 
optimal ordering only within each group. While solving (2) for 
^V, all PV-nodes and factor elements with indices corresponding 
to PV-nodes are skipped. Zeroing of reactive power mismatch at 
PV-nodes is not required. Apperantly, 50% reduction In the matrix 
storage is achieved in the absence of PV-nodes. 

PV-node Q-limit Adjustment for the FSDL method 

When the power flow solution is moderately converged, any Q- 
limit violations should be corrected. Any violated PV-node is 
converted to the PQ-type with the rective generation set at the 
limiting value. The voltage of the converted node is subsequently 
compared to the scheduled value and the node is reconverted to 
the PV-type, if any of the back-off conditions are satisfied. 

The conventional node-type switching approach is more 
suitable for the handling of the PV-node Q-limiting. The 
switching of node-types involves the insertion/delation of 
egalities into/from (2) . This can be achieved efficiently by 
partial matrix ref actorlzatlon and factor v updating methods. 
Partial Refactorization Method-1 (PRl) [13] can be used to update 
[Y"] in order to account for changes in the node status (regu- 
lated or nonregulated) . The efficiency of PRl is inversely 
proportional to the number of node status changes. Two efficient 
techniques of the node-type switching are described In the 
following. 
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(1) When network shunts are ignored (FK=0.0) or accounted by 
alternative ways, a reactivated row/ column by the removal of the 
large diagonal corresponding to a PV-node in [¥"] is identical 
to that in [Y']. It is this feature of the FSDL method exploited 
to enhance the efficiency of the node type switching using the 
following procedure. This is convenient with the simple variant 
of the FSDL method mentioned before. 

1. find the composite path of PV-nodes switched to PQ-type 

2 . find the composite path of the rest of the PV-nodes 

3. ref actor only rows/columns of [Y"] that constitute common 
path of steps 1 and 2. The factored rows/columns of [Y'] in 
the uncommon path of step 1 can be used for the partial 
ref actorization of [Y"] 

4. while solving (2), use factors from [Y'] in the uncommon 
part of the path of step 1 along with other factos of [Y"]. 

When all the PV-nodes of the system are switched to PQ-type, 
[Y'] is used for solving (2) and no partial ref actorization in 
[Y"] is required. Also when common path to be determined in step- 
3 does not exist, partial ref actorization is not required. 

(2) The compact storage scheme (1) for the FSDL method described 
above inherently provides for efficient implimentation of node- 
type switching used for PV-node Q-limit adjustment. The procedure 
constitutes the following steps. 

1. PQ-nodes in the composite path of all PV-nodes are known and 

their rows/columns are factored and stored separately from 

those of [Y'] for unadjusted solution. 
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2. Ident:lfy PQ-nodes In -the composl-te pa-th of current: PV-nodes. 

3 . Identify PQ-nodes which are not common in steps 1 and 2 . 

4. Carry out Reverse Correction Partial Matrix Ref actorization 
along the composite path of PQ-nodes of step-3 on separately 
stored factors of step-1. Use large diagonals in 
rows/ columns of PQ-nodes of step— 3 . 

If no PQ-node is identified in step-3. Reverse correction Partial 
matrix ref actorization of step-4 is not required to be performed. 
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I CLAIM : 

(a) TO have invented three Super Decoupled Loadflow methods 
(FSDL, TFDUCB and TFDLBX) both in Polar coordinate and 
Rectangular coordinates formulations, their Novel versions 
(NFSDL, NTFDLXB and NTFDLBX) , and their all possible hybrid 
cobinations. Also I claim to have invented the technique of 
producing Global corrections as given by relations (45) to 
(50) . 



(b) 



specifically the invention involves the use of the following 
seven (7) items in the methods of (a): 

(i) the definitions of the gain matrices of the six-methods 
of (a) corresponding to steps-d in algorithm-2 and 
algorithm-3, for both polar and rectangular formulations of 
the loadflow problem. Though explicit definitions, of 
unsymmetrical gain matrices of Pe-subproblem are not given, 
they are obvious for all the six methods and they are also 
claimed by this inventor as simple modifications. 

(ii) restriction of rotation angle /3 to the maximum of -36 

degrees (any other angle) in relations (3), (4), (45) and (46) 
from nonlinearity considerations, for both polar and 
rectangular coordinate formulations of the loadflow problem. 

(iii) modification of real power mismatch at PV-nodes accor- 
ding to relations (12) and (15) in FSDL, TFDLXB and TFDLBX- 
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methods formulated in both polar and rctangulr coordinates. 
This in general is the technique to modify the known vector 
to ensure the coefficient/ gain matrix symmetrical. It can 
be used in all possible computer algorithms (including all 
other approaches of the formulation of the loadflow problem) 
and real-time/ on-line/of f-line applications. 

(iv) computation of angle corrections and updating for PV- 
nodes along with the Voltage magnitude corrections while 
solving equation (24) for NFSDL, NTFDLXB and NTFDLBX 
methods. In the back-substitution part of the solution of 
(24) , I claim to have invented the technique to skip all 
the rows and columns corresponding to PV-nodes for calcula- 
tion of voltage magnitude corrections of PQ-nodes. This also 
applies when the novel methods are formulated in rectangular 
coordinates . 

(V) the Slack-start procedure for all the decoupled loadflow 
methods and in particular for the six methods of (a) above. 

(vi) similar decoupled methods as of (a) can also be used 
for solving simultaneous equations appearing in other 
areas of analysis, operation and control. 

(vii) The six invented methods of unadjusted loadflow 
solution can be used for adjusted solution by adding some 
more steps in the algorithms, state estimation, contingency 
analysis and in variety of advanced power network analysis 
and control. 
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(c) Also 1 claim t:o have invented two compact storage schemes 
for gain matrices of the FSDL method and two efficient 
procedures of the node-type switching implementations for 
the FSDL method as described in the above. 



Dated this 7 th day of november 1993. 



Signature of the Inventor S.B.Patel 
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